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ABSTRACT 

Rapid radio intra-day variability (IDV) has been discovered in the sonthern 
qnasar PKS 1257—326. Flnx density changes of np to 40% in as little as 45 min- 
ntes have been observed in this sonrce, making it, along with PKS 0405—385 and 
J1819-I-3845, one of the three most rapid IDV sonrces known. We have monitored 
the IDV in this sonrce with the Anstralia Telescope Compact Array (ATCA) at 
4.8 and 8.6 GHz over the conrse of the last year, and hnd a clear annual cycle 
in the characteristic time-scale of variability. This annnal cycle demonstrates 
nneqnivocally that interstellar scintillation is the canse of the rapid IDV at ra¬ 
dio wavelengths observed in this sonrce. We use the observed annual cycle to 
constrain the velocity of the scattering material, and the angular size of the scin¬ 
tillating component of PKS 1257—326. We observe a time delay, which also shows 
an annual cycle, between the similar variability patterns at the two frequencies. 

We suggest that this is caused by a small (~ 10 /ras) offset between the centroids 
of the 4.8 and 8.6 GHz components, and may be due to opacity effects in the 
source. The statistical properties of the observed scintillation thus enable us to 
resolve source structure on a scale of ~ 10 microarcseconds, resolution orders of 
magnitude higher than current VLBI techniques allow. General implications of 
IDV for the physical properties of sources and the turbulent ISM are discussed. 
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1. Introduction 

Some flat spectrum extragalactic radio sources display variability on time-scales shorter 
than a day. This phenomenon is known as intraday variability (IDV), and has been studied 
in detail at cm-wavelengths for two decades (Heeschen 1984; Witzel et ah 1986; see review by 
Wagner & Witzel 1995). There is now a great deal of evidence that interstellar scintillation 
(ISS) is the principal mechanism causing IDV to be observed at cm wavelengths (Jauncey 
et ah 2000; Dennett-Thorpe & de Bruyn 2001, 2002; Rickett et ah 2001, 2002; Jauncey & 
Macquart 2001). In turn, ISS can be used to probe microarcsecond (/ras) source structure and 
properties of the Galactic interstellar medium (ISM) (Macquart & Jauncey 2002; Rickett, 
Kedziora-Chudczer & Jauncey 2002). 

PKS 1257—326 is a flat-spectrum, radio and X-ray quasar at redshift = 1.256 (Perlman 
et ah 1998). Its Galactic coordinates are / = 305.2°, b = 29.9° and ecliptic coordinates 
/eel = 207.3°, feed = —24.3°. Radio IDV was discovered at 4.8 and 8.6 GHz in 2000 June 
with the Australia Telescope Gompact Array (ATGA) during blazar monitoring observations. 
These observations revealed a compact source showing strong flux density changes over the 
course of 12 hours. IDV was conhrmed in 2000 September. Here we present the results 
of observations of PKS 1257—326 carried out during a year-long monitoring program of a 
number of the known southern IDV sources (Kedziora-Ghudezer et ah 2001A; Bignall et 
ah 2002) which commenced in 2001 February. Results for the other IDV sources will be 
presented later. The observations and results are described in Section 2. In Section 3 the 
analysis of scintillation parameters and source characteristics is presented. Discussion of the 
evidence for scintillation being the principal mechanism for radio IDV in general, and of 
the connection between scintillation and intrinsic variability, is presented in Section 4. In 
Section 5 we discuss the implications of the three known rapid (intra-hour) scintillators for 
our understanding of cm-wavelength IDV in extragalactic sources. We assume go = 0.5 and 
Hq = 65 km s“^ Mpc“^ throughout. 


2. The Observations and Results 

The ATGA flux density monitoring was carried out simultaneously at 4.80 and 8.64 GHz 
during 2001 and early 2002, with observing sessions of 48 hours, approximately every six 
weeks. The target IDV sources were interspersed with a series of compact flux density 
calibrators, with observations of each source approximately every hour or less, while they 
were above 15° in elevation. This technique results in well calibrated, accurate and reliable 
flux density measurements. For PKS 1257—326, the time interval between individual mea¬ 
surements was decreased to gain good coverage of its rapid variability. Additional shorter 
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sessions concentrating on this source alone were requested and scheduled. 

PKS 1257—326 has exhibited IDV from the time it was discovered in 2000 June. Before 
this time, the only other ATCA data taken on this source were in 1995 October. We have 
inspected these data, and although they consist of only two brief (1 minute) observations, 
separated by ~ 2 hours, the two flux density measurements differ by ~ 40%; it seems the 
source was also showing rapid IDV at this time. It therefore seems likely that PKS 1257—326 
has been showing rapid IDV for at least seven years. Figure 1 shows both the 4.8 and 8.6 GHz 
ATCA flux density measurements from 2001 March 23. The flux density changes remarkably 
rapidly, from a maximum to a minimum in less than an hour, making PKS 1257—326 the 
third such rapid IDV source known along with PKS 0405—385 (Kedziora-Chudczer et al. 
1997) and J1819+3845 (Dennett-Thorpe & de Bruyn 2000A). 

Figure 1 shows the remarkably smooth, quasi-periodic nature of the variability at both 
frequencies, which we describe by a single time-scale parameter, dehned in Section 2.1 below. 
Also noticeable is the strong correlation between the variations at 4.8 and 8.6 GHz that is 
a characteristic of many IDV sources, in particular, the other extreme variables (Kedziora- 
Chudczer et al. 1997; Dennett-Thorpe & de Bruyn 2000A). 

Figure 1 demonstrates the advantage of using the ATCA for such flux density measure¬ 
ments; at these frequencies confusion is negligible, the source is mostly unresolved, and the 
array gives effectively continuous, dual-frequency, high precision flux density measurements. 
The data displayed in Figure 1 (and subsequent hgures) are averaged over 1 minute and over 
all baselines. The thermal noise for 1-minute integrations is ~ 1 mJy, neglible compared 
with the observed variations. The variation of antenna gains with elevation and time con¬ 
tributes an error which is proportional to the flux density. This is well-determined by the 
calibrator source observations. For the data published here, typical residual proportional 
errors (peak-to-peak) are up to 1% at 4.8 GHz, and up to 2% at 8.6 GHz. 

There is also an error contribution due to extended source structure. Recently we ob¬ 
tained high-resolution, high-sensitivity VLA images of the source which show it has a weak, 
arcsecond-scale jet (Bignall et ah, in preparation). The effect of this has not yet been sub¬ 
tracted from the visibilities in the data published here, however we have determined its 
contribution. The extended structure adds to the unresolved component a small amount 
of flux density, which varies over the course of a day as the (m, v) coverage of the ATGA 
varies. The contribution is different on different baselines, but the effect tends to be heavily 
smoothed out by averaging over all baselines. Nevertheless, at 4.8 GHz the extended struc¬ 
ture causes an additional, spurious “variation” of ~ 5mJy peak-to-peak at 4.8 GHz, which 
is in fact the dominant error at this frequency. At 8.6 GHz, the peak-to-peak variation due 
to extended structure is ~ 2 mJy. 
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For the purposes of estimating the time-scale and amplitude of variability from the 
light-curves, the total errors are small compared to the rapid variations, of which the rms is, 
on average, 16 mJy or 8% of the mean total flux density at 4.8 GHz, and 13 mJy or 5% of 
the mean total flux density at 8.6 GHz. Finally, the flux density scale is set through repeated 
observations of the ATGA primary calibrator PKS 1934—638 (Reynolds 1994). Thus, the 
mean flux density of PKS 1257—326 between different epochs can be compared with better 
than 2% accuracy. 


2.1. Annual Cycle in the Characteristic Time-scale 

Figure 2 presents nine of the light curves at 4.8 and 8.6 GHz made at approximately 
six week intervals over the course of a year. Immediately apparent is the dramatic change 
in the time-scale of the flux density variations at both frequencies. From February through 
May, the flux density varies rapidly with less than one hour between excursions. In June 
the variations begin a slow down which lasts through September. November sees them 
speed up again, while by 2002 January they have returned to much the same rate as 2001 
February. The time-scales estimated from the earlier discovery observations of 2000 June and 
September, and the more recent observations of 2002 February and April, are consistent with 
this annual cycle in the time-scale of variations. The annual cycle provides clear evidence 
that the observed rapid variations are due to interstellar scintillation (ISS). The change in 
the time-scale of variations through the year is due to the changing velocity of the scattering 
medium relative to the observer as the Earth revolves around the sun. 

To quantify the observed changes we dehne a characteristic time-scale of variability. 
Various authors have used different dehnitions for characteristic time-scale of IDV sources. 
These are generally based upon estimates of either the auto-correlation function (AGF) or 
the structure function of the time series, or may be estimated from the time series itself. In 
studies of scintillation of pulsars, the usual convention is to dehne the time-scale as the lag 
at which the AGF falls to 1/e of its maximum value (Gordes 1986). This dehnition of the 
“decorrelation time-scale”, tdc, was adopted by Macquart & Jauncey (2002) in their paper 
on Earth Orbit Synthesis. Dennett-Thorpe & de Bruyn (2000A) dehne the time-scale by use 
of the structure function (e.g. Simonetti, Gordes, & Heeschen 1985). Their dehnition of the 
variability time-scale corresponds to the full-width at 1/e of the AGF. Rickett et al. (1995) 
and Rickett, Kedziora-Ghudczer & Jauncey (2002) dehne the time-scale as the half-width at 
half maximum of the AGF, which we call fo.s, and which is in general slightly smaller than 

tdc- 

For extragalactic sources, in addition to the above, several other “counting” parameters 
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have been used in the literature. This is possible because the observed IDV is often quasi- 
sinusoidal, which leads to a deep minimum in the ACF, as shown in Figure 3. Kedziora- 
Chudczer et ah (1997) used the mean peak-to-peak time, while Jauncey & Macquart (2001) 
used the mean peak-to-trough or trough-to-peak time. While estimates from the ACF or 
structure function are more reliable, “counting” methods remain useful, particularly where 
the original data are unavailable. 

We calculate the ACFs using the discrete correlation method and binning into time 
lag intervals (Edelson & Krolik 1988), rather than interpolating the data. The discrete 
correlation method has the advantage that it is straight-forward to combine data from two 
(or more) consecutive days without interpolating across 12 hour gaps. The mean is first 
subtracted from each data point and the ACF is normalised to unity at zero lag. Specihcally, 
the method is as follows. For each pair of data points i,j we calculate the lag (ti — tj) and 
the function Cij = [(5'(fj) — S){S(tj) — .S)]/(t|, where S(ti) is the flux density at time U, S 
is the mean flux density for the particular dataset, and is the corresponding variance of 
S for the dataset. The elements Cij are then binned according to lag, and the average Cij 
is taken as the estimate of the ACF for each bin, provided a large number of elements he 
within that bin. Figure 3 shows ACFs calculated from the data of 2001 March 23. 

While choosing a dehnition for the time scale is a matter of convention, it is neces¬ 
sary to use a consistent method when comparing theoretical and observational results. For 
PKS 1257—326 we have calculated various time-scale estimates for each frequency and for 
each dataset, and hnd, on average, that t^c = (1.18 ± 0.05) x fo. 5 ; applies over the whole 
year and at both frequencies. For direct comparison with the results of Rickett, Kedziora- 
Chudczer & Jauncey (2002) and Rickett et ah (1995), we choose to dehne the time-scale as 
to. 5 , the 50% decay time of the ACF. 

Figures 4a and 46 show to .5 determined for each session and at both frequencies plotted 
against day of year, and show quantitatively the remarkable annual cycle so apparent in 
Figure 2 . Also included are the to .5 values calculated for the data from 2000, and the more 
recent observations of 2002 February, and April. The values for 2000 are plotted as limits. 
The light curve from June 2000 is somewhat undersampled and therefore gives an upper limit 
to to. 5 , while the 2000 September data show a decrease in flux density over 4 hours, with a 
subsequent 4.5 hour gap in the data (due to observations for another project) followed by an 
increase over the next 3 hours, so we have only a lower limit on to. 5 . Plotted in Figure 4c is the 
expected scintillation speed vs day of year, assuming that the ISM is moving with the local 
standard of rest (LSR) (Rickett et ah 2001, and private communication). Figure Ad displays 
the corresponding scintillation velocity projected onto the plane of the sky. Comparison of 
Figures 4a and 46 with Figure 4c demonstrates that the phasing of the annual cycle closely 
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matches the Earth’s projected speed with respect to the LSR. The annual cycle in to .5 is a 
result of the changing velocity of the scintillation pattern across the observer’s line-of-sight, 
and leaves no doubt as to the origin of the intra-day variability in this source. Such an 
annual cycle has been seen in two other sources; 0917-1-624 (Rickett et al. 2001; Jauncey & 
Macquart 2001), and J1819-I-3845 (Dennett-Thorpe & de Bruyn 2001), and is expected if the 
IDV is in fact a propagation effect caused by a relatively local Galactic “screen”. To reflect 
the ISS nature of the rapid variations, we hereafter refer to the characteristic time-scale, 
estimated as the 50% decay time of the ACF, as fiss- 

Scintillation is a stochastic process, and since our observations sample this process over 
a hnite duration, it is only possible to obtain an approximation to the true time-scale of 
the underlying process. This stochastic process, and our lack of understanding of it, are by 
far the dominant contributions to the uncertainty in any determination of the scintillation 
time-scale, tiss- The uncertainty decreases with the number of independent samples of the 
scintillation pattern, or “scintles”, observed. 

To estimate the error empirically, we divided the datasets from 12 well-sampled epochs, 
when tiss is short and/or the dataset extends over two consecutive days, into subsets, and 
calculate fiss for each subset. We assume that the statistical error due to the hnite number 
of samples of the stochastic process scales as 1/\/N, where N is the number of independent 
samples. We dehne N in this case as Tohs/tiss, where Tobs is the length of the observation. We 
then assume that the fractional error, (Tt/tiss, is proportional to A/tiss/T/bs- For each subset 
we calculate Et = v^[tiss(full dataset) — tiss(subset)], where the ^/2 arises because the two 
estimates are not independent. From the observed distribution of EtVN/tiss, we estimate 
the la fractional errors to be at/tiss = 0.7A/fiss/Fobs at 4.8 GHz, and at/hss = 0.9A/tiss/Fobs 
at 8.6 GHz. The resultant error bars are shown in Figures 4a and 46. 

To evaluate the reliability of the above error estimates, we assume that fiss is constant 
between day of year 50 and day of year 150 (late February through May), an assumption 
supported by the data in Figure 4, and by the almost constant value of uiss (for a medium 
moving with the LSR, as shown in Figure 4c) over this period. The observed rms scatter in 
the 8 data points from this period is 12% at 4.8 GHz and 19% at 8.6 GHz, compared with 
the mean of the errors determined for these data as above, of 13% at 4.8 GHz and 15% at 
8.6 GHz. This indicates that our error estimates, as derived above, are reasonable. 

Figure 2 also shows that the IDV at 4.8 and 8.6 GHz remained strongly correlated over 
12 months of these observations; in fact we have observed correlated variability between 
the two frequencies in every observation since the discovery of IDV in 2000 June. The 
increasing difference between the mean flux density levels reveals that the spectral index 
became increasingly inverted as the mean flux density at 8.6 GHz increased by 70%, while 
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that at 4.8 GHz increased by 35%. This is clearly seen in Figures 5a and 56, which show all 
data from well-sampled epochs at both frequencies, and Figure 5c, which shows the mean 
spectral index at each epoch. The source appears to be undergoing a pc-scale “outburst” 
of the type commonly seen in many flat-spectrum AGN (e.g., Kellermann & Pauliny-Toth 
1968). 

Figure 5d shows that despite the increase in total flux density, there is little evidence 
for a systematic increase or decrease of the rms variations at either frequency. If the increase 
in flux density were due to an increase in the flux density of the scintillating component, 
we would expect to see either a large increase in the rms variations, or if the scintillating 
component were expanding, we may observe an increase in the characteristic time-scale (in 
addition to the annual cycle). However, neither has been observed. Regression analysis 
shows that there is at best a very weak correlation in the sense that there may be a slight 
increase with time of the rms variations at 8.6 GHz. Given the low signihcance of this result, 
we are continuing to monitor the variations at both frequencies to look for any change in 
the scintillation behaviour. The connection between ISS and intrinsic variability is further 
discussed in Section 4. 

Finally, close inspection of Figure 2 reveals that at each epoch, the 8.6 GHz variability 
pattern appears to lead the 4.8 GHz light-curve. The time delay between the variations 
at the two frequencies, T( 4 , 8 ^ 8 . 6 ); was quantihed by cross-correlating the 4.8 and 8.6 GHz 
observations at each epoch. Although the overall variability patterns at both frequencies 
are quite similar, there is some extra “structure” in the 8.6 GHz light curves which appears 
to be more smoothed in the 4.8 GHz light curves. Therefore only light curves with well- 
sampled peaks and troughs give a reliable time delay measurement. As shown in Figure 6, 
the peak of the cross-correlation always has the same sign, i.e. the 8.6 GHz variations are 
leading. Although the uncertainty in each individual measurement is large, due to limited 
sampling and the fact that the light curves at each frequency are not identical, this result 
is signihcant. Interestingly, there is a clear correlation of r( 4 , 8 , 8 . 6 ) with the inverse of the 
predicted ISM velocity (Figure 4c). This gives a second annual cycle in the variability of this 
unusual source. In the next section we examine in detail the implications of these results for 
PKS 1257—326 and the ISM in our line-of-sight to this source. 


3. Microarcsecond Source Structure 

The clear presence of an annual cycle in the time-scale of variability identihes interstellar 
scintillation as the mechanism responsible for the IDV in PKS 1257—326, which in turn allows 
a determination of both source and scattering screen properties. Over the course of a year. 



the Earth’s orbital motion produces changes in both the magnitude and direction of the 
net scintillation velocity, vjss, which both potentially contribute to changes in tiss- Figure 
Ad shows the Earth’s velocity relative to the LSR over the course of a year, projected onto 
the plane transverse to the source line-of-sight. If the scintillation pattern is isotropic, only 
the change in the nett scintillation speed contributes to variations in tiss- In this case, the 
observed fiss is inversely proportional to |viss|- However, the length scale of the scintillation 
pattern may change with orientation, in which case fjss will also change as the direction 
of viss changes. An elongated scintillation pattern could be produced by elongated source 
structure and/or by anisotropic scattering in the ISM. 

Rickett, Kedziora-Chudczer & Jauncey (2002) showed that the deep first minimum, 
which they term a negative “overshoot”, in the ACEs for data on PKS 0405—385 during 
its episode of fast scintillation, indicates that the scattering is highly anisotropic. Their 
results from numerical simulation suggest that the scintillation pattern of PKS 0405—385 
has an axial ratio of 0.25 or less, where vjss is approximately aligned with the short axis 
of the irregularities. For more highly anisotropic irregularities, the depth of the ACE first 
minimum does not change signihcantly. The negative overshoot is a measure of the spectral 
purity of the fluctuations: the closer to sinusoidal the variations are, the closer to —1 is the 
depth of the first minimum in the ACE. Rickett, Kedziora-Chudczer & Jauncey (2002) argue 
that the negative overshoot is an effect of anisotropic scattering in the ISM, and cannot be 
duplicated by anisotropic source structure. 

A large number of independent samples of the scintillation pattern are required to 
determine a reliable ACE. For light curves which sample only 1 or 2 scintles, the ACE is 
poorly determined and thus a spurious negative overshoot can occur. For all of the long 
datasets on PKS 1257—326, the average depths of the first minima in the normalized ACEs 
are (—0.5 ± 0.1) and (—0.3 ± 0.1), at 4.8GHz and 8.6GHz respectively. Assuming that 
the model of Rickett, Kedziora-Chudczer & Jauncey (2002) applies to PKS 1257—326, this 
indicates that the scattering is anisotropic, with the minor axis of the scintillation pattern 
aligned within ~ 45° of the plasma velocity vector. The minima at 4.8 GHz are, on average, 
deeper than at 8.6 GHz, although there is substantial scatter (for example. Figure 3 shows 
that for the 2001 March 23 data, the hrst minimum of the ACE was in fact deeper at 8.6 
GHz.) The fact that the average size of the “negative overshoot” differs at each frequency 
may be a result of a differing source size relative to the Fresnel scale, i.e. the 8.6 GHz 
component may have a larger ratio of angular source size to angular size of the hrst Fresnel 
zone, Oq : 9^-, than the 4.8 GHz component. Hence the negative overshoot in the ACE could 
be suppressed at 8.6 GHz as source size effects become more important. 

Our data indicate that the scintillations observed at 4.8 and 8.6 GHz are in the regime 
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of weak scattering, where the medium introduces a phase change of less than 1 radian across 
a Fresnel zone (e.g. Walker 1998). Correlated variability between 4.8 and 8.6GHz is a 
pattern commonly seen in other IDV sources at high Galactic latitudes, in particular the 
other extreme IDVs. At lower frequencies, tiss is generally much longer and the variability is 
not correlated with the higher frequency variability. This is explained in terms of weak and 
strong scattering. At 4.8 and 8.6 GHz, the sources undergo weak scintillations, while lower 
frequencies (below ~ 3 GHz) are in the regime of strong (refractive) scattering (e.g. Kedziora- 
Chudczer et ah 1997; Walker 1998). Time-scales of variability for refractive scintillations 
increase with wavelength. A, as 

In the regime of weak scattering, the typical length scale of the scintillations is equal to 
the Fresnel scale, rp = a/ cL/(27rz/), where L is the distance to the scattering plasma and u is 
the observing frequency. The Fresnel scale sets an effective cut-off diameter for the angular 
size of a source undergoing weak scintillation, 9p = rp/L. A source this size or smaller 
exhibits flux density variations on a time-scale tiss = rp/viss, where viss is the speed of the 
scattering material across the line of sight. For larger source sizes, 6^ > Op, the scintillations 
decrease in amplitude as 1/Os and increase in timescale as Os (e.g. Rickett 2002). It may 
be that AGN have a sufficiently large angular diameter that the source size influences the 
observed scintillation pattern. Another effect of larger source sizes is that scintillations from 
distant scattering material are suppressed, since the angular size “cut-off”, set by the first 
Fresnel zone, scales as Therefore, there will be a bias towards hnding sources that 

scintillate behind nearby ISM turbulence. 

We can use our measurements of the annual cycle in tiss to constrain (i) the velocity of 
the scattering medium, vjsm (which Figure 4 suggests is not very different from the LSR), (ii) 
the source angular size, ^S) and (iii) the distance to the scattering screen, L. The parameters 
used to model the observed annual cycle in tiss ai’o as follows, vism may be offset from 
the LSR and this introduces two free parameters, the components of this offset in the plane 
transverse to the source line-of-sight. As shown in Figure 4d, the projected velocity of the 
sun with respect to the LSR in RA and Dec components is Vq = 16.5 km s“^, v^ = 10.7 
km s“^. The RA and Dec offsets, ^v^ and Svs, of the vjsm components shift the origin in 
Figure dd. Two more parameters are required to describe the axial ratio and orientation of 
anisotropy in the scintillation pattern. Finally there is an overall scaling factor, sq, which 
sets the scintillation length scale. We choose sq to be the minor axis of the scintillation 
pattern. We have tested the effect of varying these parameters, and hnd that the shape of 
the resulting annual cycle is quite sensitive to variations in the ISM velocity and the angle 
of the elongation for an anisotropic scintillation pattern. We hnd the best ht parameters at 
each frequency by minimising the sum of squares of the residuals. We have examined hts 
both for the case of unweighted residuals, and for the case of residuals weighted by the error 
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associated with each point. The shape of the slow-down period is critical in determining the 
ISM velocity and anisotropy in the scintillation pattern. However the errors in tiss are larger 
in this period, due to having fewer independent samples of the scintillation pattern. By using 
the unweighted residuals (i.e. ignoring the errors) we allow the £t to be better constrained 
by these points in the slow-down period. On the other hand, the overall scintillation length- 
scale is more accurately estimated using the better-sampled data from the fast period. We 
find that for models which fit the data well, using weighted or unweighted residuals makes 
little difference to the results. 


If we assume isotropic scattering and an isotropic source, the best fit to the peak in 
tiss is achieved with 6 va = —T km s“^ and 6 vs = —2 km s“^ at 4.8 GHz, where Sva and 
^v^ negative have the effect of moving the origin in Figure Ad further away from the ellipse. 
At 8.6 GHz the best fit is for 6 va = —7 km s“^ and Svs = +2 km s“^. Assuming the 
same scattering material is responsible for the scintillation at both frequencies, we take the 
average as the overall best fit: Sva = —7 km s“^ and Svs = 0. These fits, and fits for 
several other velocity offsets, are shown in Figures 7a and 7b. The scintillation length-scale. 
So, is a function of L and 63 - While we expect that source size influences the scintillation 
in PKS 1257—326, at least to the extent that scintillations from more distant scattering 
material are suppressed, the size of the scintillating component, 63 , cannot be very much 
larger than the Fresnel scale at the distance of the scattering plasma, otherwise the source 
would not show such deep modulations in flux density. The average ratio of tiss at the two 
frequencies is tiss(4.8 GHz)/tiss(8.6 GHz) = 1.3 ±0.2, which is consistent with being equal to 
0f(4.8)/0f(8-6) = a/8.6/4.8 = 1.34. Therefore the source size may in fact be smaller than, or 
equal to, the Fresnel scale at both frequencies. Alternatively, the source size may scale with 
frequency in a similar way, i.e. ^s(^) oc To consider both possibilities, we introduce a 

scaling parameter 


r 1, if 

\6*s/6*f, otherwise. 


and approximate the scintillation length scale as Sq ~ T^s ~ Re^F = Re ^/cLJ^ttu. Thus 


L ^ 0.068 



■So Y 
104 kmy 


IGHz 


pc 


( 1 ) 


Because L oc s/ uncertainties in sq lead to greater uncertainties in L, and the screen 
distance is rather weakly constrained. For the isotropic case, the allowable range of values 
found for the scintillation length-scale are sq = 6 ± 1 x 10^ km at 4.8 GHz, and sq = 5 ± 1 x 
10^ km at 8.6 GHz. Letting Rq = 1 (i.e. Sq = tf) gives an approximate upper limit on L. 
Inserting these values we find L 12±4pc at 4.8 GHz, and L ~ 15 ±6 pc at 8.6 GHz. This 
implies that the scattering occurs in a very local region of turbulence, as has been suggested 
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for the other two intra-hour scintillators, J1819-I-3845 (Dennett-Thorpe & de Bruyn 2000A) 
and PKS 0405—385 (Rickett, Kedziora-Chudczer & Jauncey 2002). If R 0 > 1, this implies 
an even closer screen. For a distance in the range 10-15 pc, r-p corresponds to an angular 
size in the range 30-37/ras at 4.8 GHz, and 22-28 /xas at 8.6 GHz. 

To obtain limits on the source brightness temperature, Tb, we also need to know the 
flux density. Sc, of the scintillating component of PKS 1257—326. Following the arguments 
of Rickett, Kedziora-Ghudczer & Jauncey (2002), Sc cannot be more than the total average 
flux density and cannot be negative, so we have St > Sc > S^ — Scam, where St is the total 
average flux density and Amin is the lowest flux density observed. Since we know At varies 
with time, we compute At —Amin for each epoch, and hnd that the maximum value is close to 
50 mJy at both frequencies. Also, At at the start of our monitoring was ~ 200 mJy at both 
frequencies. In fact we can place a further constraint on Ac at 4.8 GHz, since a short-baseline 
VLBI observation of the source, including the ATGA and an unresolved calibrator to set the 
flux density scale, showed that ~ 25% of the total average flux of PKS 1257—326 (on 2001 
March 23) is missing on an angular scale of ~ 0.1" at 4.8 GHz. Thus we assume 50mJy 
< Ac < 150 mJy. 

Assuming a circularly symmetric source brightness distribution with FWHM corre¬ 
sponding to the Fresnel scale size, an approximate lower limit on the rest-frame brightness 
temperature, Tb = 4 x 10^^ K, is obtained for Ac = 50mJy and 9s = 40/ias at an observed 
frequency of 4.8 GHz. This requires relativistic beaming with a Doppler factor of "D > 4 
so as not to violate the Inverse Gompton (IG) limited brightness temperature of 10^^ K for 
synchrotron radiation. If in the rest frame the source is at an equipartition brightness tem¬ 
perature of ~ 2 X 10^^ K, as defined by Readhead (1994), then T ~ 20 is required. While 
high, these implied Doppler factors are not inconsistent with apparent superluminal speeds 
found in VLBI surveys (Marscher et al. 2000; Kellermann et al. 2000). Taking Ac = 150 mJy 
and 6*s = 30/ias gives a somewhat more extreme brightness temperature of Tb = 2 x 10^3 K. 

Figures 7c and 7d show hts obtained allowing an anisotropic scintillation pattern with 
no velocity offset from the LSR. As discussed above, there is independent evidence for 
anisotropic scattering, based on the depth of the hrst minimum in the AGFs. Finally, 
Figures 7e and 7/ show hts obtained for an anisotropic scintillation pattern and allowing 
a screen velocity offset from the LSR. The width of the tiss peak, allowing all parameters 
to vary, is best htted for both frequencies with Jvq, = —5 km s“^ and Svs = 0, and when 
the minor axis of the scintillation pattern lies at an angle of approximately 25° ± 10° to 
visM- The ISM velocity offset is shown by the arrow in Figure Ad, and the orientation of 
the htted scintillation pattern is shown by the small ellipse in Figure Ad. The best hts are 
found for axial ratios < 0.5. Decreasing the axial ratio below ~ 0.25 has little ehect on 
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the shape of the tiss annual cycle, because vjss never cuts directly across the long axis of 
the htted scintillation pattern. This result is consistent with what is expected from depth 
of the hrst minimum in the ACFs, which in the case of PKS 0405—385 was htted with a 
highly anisotropic medium (axial ratio ~ 0.25; Rickett, Kedziora-Chudczer & Jauncey 2002). 
Those sources for which viss at its maximum cuts across the short axis of the scintillation 
pattern are those most likely to scintillate rapidly. 

For the best hts allowing all parameters to vary, the minor axes of the scintillation 
pattern scale are approximately Sq = 4.2±0.5 x 10^ km at 4.8 GHz, and Sq = 3.5±0.5 x 10^ km 
at 8.6 GHz. In the case of anisotropic scattering, the typical length scale for the minor axis 
is the Fresnel scale reduced by a factor of 72 (B.J. Rickett 2002, priv. comm.). We 
estimate L by setting Rg = 1 and replacing sq with r-p ^ 72so in Equation 1. This gives 
values for the screen distance of L ~ 11.5±3pc at 4.8 GHz, and L k, 14.4±4pc at 8.6 GHz; 
very similar values to those obtained for the isotropic case. Hence the implied brightness 
temperatures are also very similar to those found for the case of an isotropic scintillation 
pattern. However, in the case of anisotropic scattering, it may be that Tb could be reduced 
further if the source was coincidentally elongated in the same direction as the long axis of 
the scattering. 

For rapid scintillators such as PKS 1257—326, it is possible to measure the IDV pattern 
time delay between two widely separated telescopes (Jauncey et al. 2000; Dennett-Thorpe 
& de Bruyn 2002). Recently we observed PKS 1257—326 simultaneously using the ATGA 
and the Very Large Array (VLA), and observed a clear time delay between the IDV pattern 
arrival times at each telescope (Bignall et ah, in preparation). By measuring this delay 2 or 
3 times over the course of a year, we will obtain an accurate, independent estimate of the 
scintillation parameters. 


3.1. Annual cycle in the time offset betw^een 4.8 and 8.6 GHz 

We now turn our attention to the annual cycle in r( 4 . 8 , 8 . 6 ); the offset between the IDV 
patterns observed at 4.8 and 8.6 GHz. We suggest that this may be due to an offset between 
the central components of the source at each frequency, as might be expected if the source 
were jet-like on a /xas scale, and optically thick between 5 and 8 GHz. The line in Figure 6 
shows an expected annual cycle for an offset between the two components, with projected 
displacement vector D{a,6) = (1.5, 0.9) x 10^ km in components of RA and Dec, using the 
best-ht ISM velocity derived from the tiss annual cycle above. Given the large (but difficult 
to estimate) uncertainties in the measurement of T( 4 . 8 , 8 . 6 )) fhis model seems a reasonable ht 
to the data. The fact that T( 4 , 8 , 8 . 6 ) always has the same sign means that the ISM always 
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crosses the 8.6 GHz “component” first, which severely constrains the direction of the offset 
of the two components, and also implies that the origin in Figure Ad cannot be inside the 
ellipse; i.e. viss changes direction by < 180° over the year. This is consistent with what is 
expected from the annual cycle in tiss- 

For a screen distance of 10 pc, the fitted displacement vector corresponds to an offset 
of 12/ias. This order of resolution is currently unachievable by any other technique, and 
is comparable to that expected for proposed future space-based X-ray (Cash et ah 2000) 
and optical (Unwin et al. 2002) interferometry missions. An offset of this order between 
different frequency components may also have implications for future astrometry and geodesy 
programs. An angular scale of 12/ias at z = 1.256 corresponds to a linear scale of 0.08 pc. 

Another possibility which we have not ruled out is that the frequency offset is due to 
a refractive effect in the ISM rather than the source itself. Theoretical investigation of such 
effects is ongoing and will be presented in a later paper. 


4. Discussion 

The evidence for interstellar scintillation as the principal cause of IDV at cm wavelengths 
is now compelling. Three sources have shown a clear annual cycle in the changes in their 
characteristic time-scales, J1819-I-3845, 0917-1-624, and now PKS 1257—326. In addition, 
a time-delay has been measured in the arrival times of the IDV pattern at widely spaced 
radio telescopes for PKS 0405—385 (Jauncey et al. 2000), J1819-I-3845 (Dennett-Thorpe & 
de Bruyn 2002), and recently PKS 1257—326 (Bignall et ah, in preparation). 

At this point, it is important to re-examine the evidence presented for an intrinsic origin 
for cm-wavelength IDV. The strongest case reported is the correlated radio-optical changes 
found in 0716-1-714 (Quirrenbach et al. 1991; Wagner et al. 1996). In 1990 February, optical 
variations were found that were strongly correlated with changes in the radio spectral index 
between 6 and 3.6 cm wavelengths, in the sense that the optical flux increased as the radio 
spectrum became more inverted. If this correlation were to extend across the spectrum from 
radio to optical wavelengths, as would be the case if the changes were intrinsic, then such a 
spectral inversion should lead directly to even stronger flux density changes at wavelengths 
shorter than 3.6 cm. However, such changes are not evident in the 2 cm data for this source 
from the same epoch (Quirrenbach et al. 2000). This raises some doubt as to the value of 
the evidence for the intrinsic nature of the variability seen in this source in 1990 February. 
It is notable that 0716-1-714 continues to show dramatic optical variability on a time-scale of 
days or less (Nesci, Massaro & Montagni 2002). 



Alternatively, the low level of IDV in 0716+714 at both 20 and 2 cm is consistent with 
a scintillation origin, as is seen in the variations in amplitude of IDV with frequency in 
both PKS 0405—385 (Kedziora-Chudczer et ah 1997) and J1819+3845 (Dennett-Thorpe & 
de Bruyn 2000A). Moreover, the apparent difference in phase between the 5 and 8 GHz 
data resembles that apparent in PKS 1257—326, as well as in J1819+3845 (see hgure 2 of 
Dennett-Thorpe & de Bruyn 2000A), further strengthening the case for a scintillation origin 
of the radio IDV in 0716+714. Given its proximity on the sky to 0917+624, it is important 
to search for an annual cycle in the variability time-scale of this source, similar to that seen 
in 0917+624 (Rickett et al. 2001; Jauncey & Macquart 2001). 

The case for interstellar scintillation at cm-wavelengths is so strong that we suggest that 
this phenomenon be referred to as the physical phenomenon of interstellar scintillation, rather 
than the observational phenomenon of IDV. This is not merely a change in nomenclature, but 
has value in focusing on the physics; while IDV is scintillation, scintillation is not necessarily 
IDV. For example, in the strong (refractive) scattering regime, scintillating sources show a 
characteristic time-scale which increases with decreasing frequency (oc that is easily 

an order of magnitude or more longer than intra-day. 

What is the connection, then, between scintillation and intrinsic variability at cm wave¬ 
lengths? Variability of compact, flat-spectrum radio sources has been studied in detail for 
four decades. These sources often show large “outbursts” on time-scales of months to years 
at cm wavelengths, which have been successfully modelled as intrinsic variations in relativis¬ 
tic jets, and linked with apparent superluminal motions of components observed with VLBI 
(e.g. Kellermann & Pauliny-Toth 1981). While refractive scintillations were found to be im¬ 
portant for low frequency variability (frequencies < 1 GHz) (Hunstead 1972; Rickett, Goles, 
& Bourgois 1984), for many years the prevailing paradigm was that cm wavelength variabil¬ 
ity was all intrinsic. It was only with the discovery of rapid IDV (e.g. Kedziora-Ghudczer 
et al. 1997) implying excessively high brightness temperatures if the variability was intrinsic 
(10^^ K in the case of PKS 0405—385), that the importance of ISS at cm wavelengths has 
been extensively re-examined. In particular, the largest intensity fluctuations due to ISS are 
observed at frequencies close to the transition between weak and strong scattering, typically 
~ 3 — 8 GHz along lines of sight out of the Galactic plane (Walker 1998). 

Our flux density monitoring of PKS 1257—326 suggests that it is also undergoing the 
type of intrinsic, pc-scale outburst commonly seen in flat-spectrum sources, yet there has 
been no signihcant change in the scintillating flux density. This is despite the optically thick 
spectrum of the slowly varying source, and suggests that the /ras core does not lie behind 
what is presumably an expanding pc-scale jet, as might be expected if the scintillating 
component were the core at the base of the pc-scale jet. In other words, the outburst is 
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probably occurring in a separate emitting region, further out from the core in a component 
too large to scintillate. High resolution VLBI observations could provide information on 
activity in the source on a scale of several pc. 

It is interesting to compare the brightness temperature inferred from the long-term 
intrinsic variability, with that obtained from the ISS model. For PKS 1257—326, we observed 
a steady increase in the mean flux density at 8.6 GHz of AS = 0.14 Jy over a period of 440 
days, while over the same time the mean 4.8 GHz flux density increased by AS = 0.08 Jy. 
The simple approach of taking these values of AS to be the flux density of the slowly varying 
component at each frequency, and observed variability time-scale tint = 440 days, gives a 
variability brightness temperature (in the source proper frame), T™'’, of 3 x 10^^ K from the 
8.6 GHz long-term variability, and T™'’ = 6 x 10^^ K from the 4.8 GHz long-term variability 
(Lahteenmaki & Valtaoja 1999). For relativistically beamed radiation with Doppler factor 
V, the inferred variability brightness temperature T™ is related to the source rest-frame 
brightness temperature T^, by T™ = V^TS Where Tb is calculated from an angular size 
which has been measured directly, or estimated from an ISS model as in Section 3 above, 
then Tb = 'DTS Taking this dependence into account, the intrinsic brightness temperature 
inferred from the long-term variability is smaller than that calculated from the modelled 
angular size of the scintillating component. This is consistent with the long-term variability 
occurring in a larger region, separate from the scintillating component, and the lack of any 
detectable change in the scintillating flux density. 

In some sources, variability due to ISS is episodic. For example, the IDV in PKS 0405—385 
has been convincingly shown to be due to ISS (Kedziora-Ghudczer et al. 1997; Jauncey et 
ah 2000; Rickett, Kedziora-Ghudczer & Jauncey 2002), however this source has shown only 
short-lived episodes of fast scintillation, each lasting several months (Kedziora-Ghudczer et 
al. 2001B). PKS 0405—385 also exhibits slow flux density evolution, and Kedziora-Ghudczer 
et al. (1997) proposed a possible interpretation of the episodic ISS, relating the long-term 
behaviour to the ejection and subsequent expansion of a component small enough to scin¬ 
tillate. However, this interpretation did not agree so well with the second episode of ISS 
observed in 1998, and it remains unclear whether the transient nature of the scintillation 
observed in PKS 0405—385 is due to intrinsic changes in the source occurring on time-scales 
of months, which cause the appearance and disappearance of a scintillating component, or 
to changes in the properties of the scattering medium. 

A similar problem exists for 0917-1-624. Recent monitoring data on this source were 
presented at the AGN Variability Workshop in Sydney by Fuhrmann et al. (2002), showing 
that the rapid, large-amplitude scintillation of this source has recently ceased after 15 years 
of monitoring. More recent VLA data collected during the MASIV Survey (Lovell et al. 
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in preparation), and data from the Goldstone-Apple Valley Radio Telescope program (M. 
Klein, priv. comm.) confirm that the rapid scintillations in 0917+624 have so far not 
returned. 

While ISS and observations of annual cycles are useful for probing pas-scale source 
structure, there is clearly much remaining to be understood about ISS of extragalactic radio 
sources, and its connection to the intrinsic variability which generally occurs on longer time- 
scales. 


5. Significance of the Fast Scintillators 

The three fast scintillators PKS 0405—385, J1819+3845, and PKS 1257—326, are valu¬ 
able because their rapid scintillation time-scale allows one to characterize the statistics of the 
underlying scintillation process over a typical 12-hour observing session; one observes more 
stastically independent “scintles” per unit time compared with slower scintillators. This is 
important because the dominant source of error in modelling the structure of a scintillating 
source is in the stochastic properties of the scintillation process itself. Moreover, in their 
periods of rapid variability it proved possible to make the pattern time delay measurements, 
crucial in establishing the ISS origin (Jauncey et al. 2000; Dennett-Thorpe & de Bruyn 2002), 
since this is not possible for the more usual, slower scintillating sources with characteristic 
time scales of order a day. These three rapidly scintillating sources proved invaluable in mak¬ 
ing the transition from IDV to ISS. Viewed from an IDV perspective, their properties appear 
extreme. From a scintillation perspective, however, are these sources really so extreme? 

The answer lies in the reasons for their rapid scintillation. The characteristic time-scale 
for scintillation is fiss = -siss/'^iss, where siss is the scintillation length scale and viss is the 
bulk velocity of the scattering plasma in the observer’s frame of rest. If an annual cycle is 
observed in tiss, it follows that the scattering plasma velocity must be close to the 30 km 
s“^ of the Earth’s orbital speed. As discussed in Section 3, for an observed frequency u in 
the weak scattering regime, siss is related to the Fresnel scale, rp, which scales with screen 
distance L as rp oc \/L. siss also depends on the source angular size, 6 s, in the case where this 
exceeds the angular Fresnel scale, such that siss ~ LOs- For most extragalactic sources, the 
relation of time scale to screen distance lies between a square root and a linear dependence 
(Rickett, Kedziora-Chudczer & Jauncey 2002). It follows that the fast scintillators, those 
with small fiss, will be those scintillating sources which are to be found behind very nearby 
scattering screens. The evidence in support of this is the screen distances of < 30 pc found 
for J1819+3845 (Dennett-Thorpe & de Bruyn 2000A), PKS 0405—385 (Rickett, Kedziora- 
Chudczer & Jauncey 2002), and PKS 1257—326. 
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From the scintillation perspective, these fast variables are extreme in the sense that they 
are to be found behind nearby scattering material, not in that there is any inherently extreme 
property in the sources themselves. There is no evidence for anything intrinsically unusual 
about PKS 1257—326 which sets it apart from other flat-spectrum, radio-loud quasars. It 
has optical magnitude B = 18.7 and has been detected at X-ray energies, emitting 2 x 10“^^ 
erg cm“^ s“^ between 0.1 and 2.0 keV (Perlman et ah 1998). In fact, its radio properties are 
less extreme than might be expected for the more common, slower scintillators. Since the 
angular size required for a source to scintillate in the ISM is set by the angular size of the 
hrst Fresnel zone, those sources which show large variations with longer characteristic time 
scales are more likely to be found behind scattering screens at much larger distances, and 
hence with smaller angular sizes and correspondingly higher brightness temperatures. 

There are also implications for the properties of the interstellar medium. The small 
number of fast scintillators suggests that there are few regions of the Galaxy where such 
nearby screens may be found, a point supported by the observation that the screen towards 
J1819-I-3845 does not move at the local standard of rest (Dennett-Thorpe & de Bruyn 2002). 
Moreover, as noted above, the presence of a nearby screen relaxes the constraints on the 
angular size necessary for a source to scintillate. Thus there will likely be more sources with 
components with angular sizes of order 20 pas, which can scintillate through a screen at 
L < 30 pc, than there will be sources with 5-10 pas angular sized components to scintillate 
through a screen at, say, 500 pc. In this regard, it is noticeable that the two long-lived fast 
scintillators, J1819-I-3845 and PKS 1257—326, are the two weakest IDV sources reported. 
If the above is correct, then a deep, unbiased scintillation survey would be expected to 
conhrm the low numbers of fast scintillators. Alternatively, the discovery of more such rapid 
scintillators can reveal the presence of such nearby scattering material in the Galaxy. A 
survey is currently underway, using the VLA to look for scintillation-induced variability over 
the whole Northern sky, in both strong (> 0.5 Jy) and weak (~ 0.1 Jy) flat-spectrum sources 
(Lovell et ah in preparation). 


6. Summary 

Here we summarize the main results from our observations of PKS 1257—326. 

• PKS 1257—326 is a long-lived, rapid IDV radio source with similar properties to the 
other rapid IDVs, in particular PKS 0405—385 (Kedziora-Ghudczer et al. 1997; Rickett, 
Kedziora-Ghudczer & Jauncey 2002) and J1819-I-3845 (Dennett-Thorpe & de Bruyn 
2000A). 



• Monitoring over more than a year reveals a clear annual cycle in the characteristic 
timescale of variability. 

• This annual cycle provides unequivocal evidence for a scintillation origin of the IDV 
observed in this source. 

• Modelling the annual cycle provides strong evidence for highly anisotropic scattering 
in the ISM. 

• The very rapid variability is most likely due to scintillation in a nearby screen, only 
~ 10 to 15 pc from the sun. 

• Assuming a source size equal to the Fresnel scale at the distance of the screen, the 
implied brightness temperature is then 4 x lO^'K < Tb < 2 X 10^^ K, which is high, 
but not entirely unacceptable for relativistically beamed radiation from a jet. 

• A time offset is observed between the scintillation patterns at each frequency, which 
itself exhibits an annual cycle. This we tentatively model as an offset due to opacity 
effects in the source on a scale of ~ lO/ras, or ~ 0.1 pc at the source. 

• There are much slower flux density changes in the pc-scale structure of the source, but 
these do not appear to be mirrored in any obvious changes in the microarcsecond sized 
component responsible for the IDV. 

• We plan to continue to monitor PKS 1257—326, to better constrain the time-scale 
during the slow-down period, and to follow pc-scale changes in the source with high- 
resolution VLBI observations. 

• For sources scintillating in the weak scattering regime which show large variations 
on timescales of the order of a day, the scattering screens may be at much larger 
distances (several hundred pc). This requires smaller source sizes, so in fact these 
slower scintillators may have higher implied brightness temperatures than the very 
rapid scintillators. 

• Observations of scintillating sources provide a method of achieving /ras resolution at 
cm-wavelengths with modest radio telescopes, resolution comparable with that pro¬ 
posed for future space-based optical and X-ray interferometers. 


The Australia Telescope Compact Array is part of the Australia Telescope which is 
funded by the Commonwealth of Australia for operation as a National Facility managed 
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Fig. 2.— Light curves over the course of a year for PKS 1257—326. All plots are on the 
same scale, shown on bottom left plot. Large symbols represent 4.8 GHz data; small symbols 
(colored red in the electronic edition) represent 8.6 GHz data. 
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Fig. 4.— Characteristic timescale to .5 (defined as the HWHM of the ACF) vs day of year at 
(a) 4.8GHz and (b) 8.6GHz. Values for 2000 are limits (see text for details), (c) Expected 
scintillation speed vjss vs day of year for a scattering medium moving with the LSR. (d) The 
corresponding scintillation velocity viss projected onto the plane transverse to the source line- 
of-sight, in components of RA (va) and Dec {vs). Arrow shows offset velocity of {6va, Svs) = 
(—5,0)kms“^. Small ellipse shows a contour of the best fit scintillation pattern. See text 
for details. 






Fig. 5.— Flux density at (a) 4.8 and (b) 8.6 GHz, from all well-sampled epochs, plotted 
with 1-minute averaging, (c) Mean spectral index from each epoch, (d) rms variation for 
each epoch at 4.8 GHz (open circles) and 8.6 GHz (closed circles). 












Fig. 7.— Various fits to the annual cycle in tiss- (a) and (b) are for an isotropic source 
and isotropic scattering, (c) and (d) are for a medium moving with the local standard of 
rest (LSR), allowing an anisotropic scintillation pattern, (e) and (f) are for an ISM velocity 
offset from the LSR and an anisotropic pattern. Model parameters are shown on the plots. 
Ratio is the axial ratio of the anisotropy. 6va and Svs are in units of km s“^, Sq is in units 
of 10"^ km. 6 is the angle between vism and the minor axis of the scintillation pattern. See 
text for further explanation. 



























